Glueball masses from an infrared moment problem and nonperturbative Landau gauge 
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We set up an infrared-based moment problem to obtain estimates of the masses of the scalar, pseudoscalar 
and tensor glueballs in Euclidean Yang-Mills theories using the Refined Gribov-Zwanziger (RGZ) version of 
the Landau gauge, which takes into account nonperturbative physics related to gauge copies. Employing lattice 
input for the mass scales of the RGZ gluon propagator, the lowest order moment problem approximation gives 
the values m 0++ ~ 1.96 GeV, m 2 ++ ~ 2.04 GeV and m -+ ~ 2.19 GeV in the SU(3) case, all within a 20% 
range of the corresponding lattice values. We also recover the mass hierarchy mQ ++ < ra2++ < m () h- 
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There is little doubt that Quantum Chromodynamics (QCD) 
displays confinement: the elementary excitations, viz. the 
quarks and gluons, which carry color charge, are not observ- 
able. They appear in the physical spectrum only in color- 
less combinations. Although confinement has not yet been 
proven, the experimental and numerical lattice evidence is 
overwhelming. Next to the hadronic sector, the spectrum of 
QCD is expected to exhibit glueballs, colorless purely gluonic 
states. Their observation is, however, cumbersome, due to the 
mixing with mesonic states with identical quantum numbers. 
This complication can be circumvented by first studying pure 
SU(A0 Yang-Mills (YM) theory, i.e. QCD without quarks. In 
this setting, there is ample lattice evidence for the existence 
of glueballs, see |T| for an overview on the glueball related 
literature. 

In the continuum approach to QCD, one usually needs to 
choose a gauge. It well known that a local and Lorentz covari- 
ant gauge, which singles out only one representant per gauge 
orbit, does not exist 0. Gribov was the first to realize that the 
elimination of gauge copies is a non-trivial issue [3]. Take, 
for example, the Landau gauge d^A^ = 0, and consider an 



(infinitesimal) gauge transformation, A'^ = A 



Then 



dfjA'^ = if d^DfjG) = 0. Thus, Landau gauge copies certainly 



occur when the Faddeev-Popov operator M" b 

->ofc _ 3 %ab gf abc > 



D" b = dfi ab — gf abc A c has zero modes. To exclude this class 
of gauge copies, it was proposed to restrict the functional 
integration to the region D. = {A^d^ = 0,M ab > 0}. Gri- 
bov, and later on Zwanziger, were able to translate this restric- 
tion into the framework of local quantum field theory J3j 0] . 
This leads to an improved Faddeev-Popov action, known as 
the GZ action, 



Sgz 



d 4 x 



-gf^d^D^fc^ + y 4 d{N 2 - I) 



(1) 
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with {<pjf,?jf} and {cof,co^} a set of bosonic, 
resp. fermionic fields. The parameter y, carrying the di- 
mension of a mass, is not free, but corresponds to the 



nonvanishing solution of the gap equation 



dy 



0, known 



as the horizon condition [4|, where E vac is the vacuum energy. 
The presence of the parameter y ~ Aqcd and the horizon con- 
dition have a deep influence on the IR behavior of the gluon 
and ghost propagator, a result which cannot be achieved with 
perturbation theory derived from the Faddeev-Popov action. 
In lEl El, some of us showed that additional effects, under the 
form of d = 2 condensates, have to be taken into account, 
giving rise to the so called Refined Gribov-Zwanziger action. 
The resulting gluon and ghost propagators are in agreement 
with lattice data at large volume [7). Here, we suffice to 
mention that a major role is played by dynamically generated 
mass terms which can be accommodated by the following 
version of the GZ action J5] |6), 



Sgz -^Sgz+I d 4 * ( ^A^; - M 2 <s?f cpf 



(2) 



which accounts for a d = 2 condensate in the gluon sector [8] 
and one in the {<p<f ,(af,W^ b } sector ||5]|6]]. The tree level 
gluon propagator (Affi) k = h ab (5^ v - k^/k 2 ) D{k 2 ) is 



D{k 2 



fc 4 + (m 2 +M 2 )fc 2 + A, 4 



X 



2g 2 Ny 4 + m 2 M 2 . 



(3) 



In (9), it was tested inhowfar Q could reproduce the lattice 
data quantitatively. An accurate fit was possible up to k ~ 
1.5 GeV, leading to the continuum extrapolated values 



i 0.337 GeV 2 ,M 2 « 2.15 GeV 2 ,V ~ 0.26 GeV 4 



As we did not consider loop corrections, it is no surprise that 
the perturb ative/UV logarithmic tail is not well-reproduced 
|9l . It is believed that relevant nonperturbative contributions 
come from momentum space gluon configuration A(k) for 
k < 1 GeV, see e.g. [10|. This suggests that expression ^ 
can capture a large amount of the nonperturbative/IR content 
of the gluon propagator. It is easily seen that expression ([3]) 
has 2 cc poles, given in GeV units by 

-k 2 = m\ = f? ± iV2Q 2 = 0. 1685 ± 0.48 12/ . 
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We introduced the parameters /j 2 and 2 for later convenience. 
These cc masses imply a violation of positivity [12] in the 
gluon propagator, again in compliance with numerical data 
ifTTl . This can be seen as a manifestation of confinement: 
gluons cannot propagate as free physical particles. We re- 
mind here that, in order to have a physical meaning, an Eu- 
clidean two-point correlation function F(k 2 ) should exhibit 
the Kallen-Lehmann spectral representation lfT2ll 



F(k 2 ) = 



TO 



p(0 

t + k 2 



dt , 



(4) 



where p(f ) has to be positive for t > To, as it corresponds to a 
cross section by the optical theorem. The spectral density p(f) 
contains a large amount of physical information. For example, 
a pole at k 2 = —m 2 in F{k 2 ), i.e. a physical particle with mass 
m*, corresponds to a delta function, 8(f — m 2 ), present in p(r). 
Usually, as YM is asymptotically free, one computes large 
momentum contributions to F(k 2 ), supplemented by nonper- 
turbative OPE contributions, which are employed to estimate 
physical quantities for which a truly nonperturbative analytic 
evaluation is unknown. For instance, a single resonance plus 
a continuum parametrization, p(f ) = aS(f — m 2 ) + bQ(t — to), 
enables one to obtain estimates of the spectrum, being at the 
basis of what is referred to as sum rules, moments problem, 
etc., see HI. 

In the current GZ context, a different approach seems to be 
more appropriate. We recall that the propagator |3]l is ideally 
suited to describe the IR sector of YM theories in the Lan- 
dau gauge. We shall first describe our set up, and afterwards 
explain why we feel that it is the most appropriate way to pro- 
ceed. Setting t = l/s in Q yields 



F(k z 



!Ao p(l/j) 



1 

l+sk 2 



ds = 



l+sk 2 



ds, 



which is ideally suited to be Taylor expanded at small k 

F(fe 2 )=£v„(-ir(fe 2 r ) 



where we defined the moments 



v„ = / s"a(s)ds. 
Jo 



(5) 



(6) 



Let us also introduce 



Z \ Z J Jo z — s 

which provides us with a more standard formulation for the 
(reduced) Hausdorff moment problem with finite boundaries 
04). In general, the (reduced) moment problem amounts to 
search for a spectral function a(s) > capable of generating 
a (finite) set of numbers v„, according to d6l. This problem 
has been solved in terms of Pade appro ximants [14|, which 
are rational approximations to a given function, which is here 
a glueball propagator. The poles of this rational approxima- 
tion can be interpreted as mass estimates for the unknown full 
propagator, of which we only know the approximation F(k 2 ). 



Usually, the Pade approximant is a much better approxima- 
tion to a function than its Taylor series. An expansion at small 
k 2 > corresponds to an expansion in ~ near z ~ — °°. By 

power counting, we see that f(z) ~ \ for z ~ — °°, thence we 
shall need a Pade approximation of the type [N,N—1] lfT4l . 



(7) 



where N refers to the order of the polynomial. By match- 
ing l.h.s. and r.h.s. of |7]) up to the designated order, the co- 
efficients of the polynomials Pjv-i (z) an d Qn{z) can be com- 
pletely determined in terms of the moments v„. The will 
be orthogonal polynomials over [0,£o] with weight a(s) lfT4l . 
Consequently, their poles z* will all be real, different and lying 
in the interval ]0,£o[. Concerning the mass estimates, these 
will be then necessarily real and larger than 1 /Eq = Tq, with 
To the threshold appearing in (ffl. An important role is played 
by the positivity of a(s) or, more precisely, by the require- 
ment that the moments v„ correspond to a positive weight 
g(s). This can also be characterized completely in terms of 
the Pade approximants, see lfT4l . 

Before proceeding, we spend a few words on subtractions. 
The UV infinities of a Green function like F(k 2 ) will be re- 
flected in the large f-behaviour of p(f), making the integral 
Q to blow up. These infinities can be dealt with by turning 
to a subtracted spectral representation. By a sufficient number 
of derivatives w.r.t. k 2 , let us say r, the integral Q can be reg- 

ulated. By integrating back each time from T to k , we 

get the following finite subtracted representation 



F sub (k z ) = (-l) r {k z -T)' 



1 



s r a(s) 
(l+sT) r l+sk 2 

S v ' 



ds. 



T is the subtraction scale, which can be chosen at will. As 
physical results should not depend on T, we shall use the lore 
of minimal sensitivity ifTSIl . Notice that acting with j^zy is 
sufficient to kill the divergent moments Vo, . . . , V r -i in the se- 
ries ((5J, the residual moments are well-defined. It is clear 
that any reasonable scheme capable of discussing the gener- 
ation of poles in the subtracted F sub (k 2 ), shall do the same 
for F(k ), since the difference between both expressions is 
a polynomial in k 2 . We shall thus consider the IR moment 
problem associated to the v' n , which we define using F(k 2 ) = 

/ r ° ^jjds, so that we can introduce f(z) = (— ~) with 



Vj, = J*W(j)«fe<oo. 



We still need to specify exactly which moment problem we in- 
tend to solve. As mentioned, the propagator Q gives a good 
description of the IR region. Evidently, a glueball operator 
constructed from such gluon propagator cannot be a true ap- 
proximation over the whole momentum region. Nevertheless, 
we might expect that it can be trusted in the IR. Indeed, at tree 
level, we can always write 



7 sub 



(k 2 ) 



{-\y(k 2 -Ty 



dq 



/(*,?) 



(8) 
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where / assembles information on the gluon propagator, 
tensor structures, etc. Since the r.h.s. of |8]l is supposed 
to be finite, so must be the l.h.s., as such we know that 
\im q ^ m f(k,q) < °° based on power counting. For small k 2 , 
the value of the integral (|8]l should be dominated by low 
q 2 , i.e. by low momentum information of the gluon propa- 
gators, since we can expand the function f(k,q) — f(Q,q) + 

^{0,q)k 2 + Generically, acting with g|j on f(k,q) can 

only bring down more powers of q 2 in the integrand, thereby 
further diminishing the effect of the UV. Similar arguments 
hold for higher loop integrals, meaning that we should only 
keep IR information in our glueball propagator. In general, 
one could also construct an UV moment problem by expand- 
ing Q for large k 2 , but in the light of the previous remark, this 
appears paradoxical, as the input of a Gribov-like propagator 
is clearly of an IR nature. We shall thus only keep the low- 
est IR moments v' ; . The more moments we keep, the "less" 
IR and hence trustworthy the data becomes. In practice, we 
shall only keep the 2 first moments Vq and v',, which usually 
already give a very good approximation to the function one 
starts with. We have now to specify which RGZ operators we 
identify with glueballs. In general, given a glueball with quan- 
tum numbers J , one needs a composite operator with ex- 
actly those numbers. In the YM context, this is done by look- 
ing at the appropriate classically gauge invariant operator. For 
the sake of presentation, let us discuss here briefly the scalar 
glueball, which ought to be related to F 2 V . We shall work 
at lowest order, thus we can look at its Abelian content f 2 v . 
Upon inspection of the action ([TJ, it is clear that the gluon and 
the new fields are partially mixed up. We already mentioned 
the occurrence of the 2 cc masses in the corresponding gluon 
propagator. It is useful to introduce the so-called /-particles 
which enable us to diagonalize the kinetic part of the action 
(|2). Using lineai- combinations of A^ and {<p^,<p^} ED, one 
can diagonalize the tree level piece of Q 



/ 



d 4 x 



-a» (-a 2 + m 2 + ) x; + -r,; (-a 2 + m 2 _ ) ^ + 



rest 



Clearly, the cc fields and r|", have the cc masses m\. 
We can introduce the field strengths A,^ v = d^k^ — and 

ri", v = - d ^K' and verif y that fly = KvK* + rest - 

Details about the foregoing construction will be reported else- 
where. The first part A-^ V T|^ V will generate a physically mean- 
ingful piece in the corresponding correlator, i.e. with a Kallen- 
Lehmann representation, while the second part is responsible 
for an unphysical piece 11161 . For the remainder of this work, 
we shall simply focus on the physical part of the operator. We 
repeat that this approach is currently only an assumption. In 
ongoing work, we are investigating whether this kind of oper- 
ator can be consistently extended to the interacting quantum 
level and whether it can be constructed in accordance with the 
softly broken BRST [6| or with the modified BRST transfor- 
mation of 02). The previous assumption about what a physi- 
cal glueball operator would look like will be justified here by 
the numerical estimates for the masses which we shall be able 
to work out. 

Let us investigate now explicitly the scalar, pseudoscalar and 



tensor glueball. For the scalar we can employ the physical part 
of F 2 V , while for the pseudoscalar we look at ie^/^vFcp. 
For the tensor, we would use the standard operator ifPIl 



a/j 1 av 



(9) 



which is the classical energy-momentum tensor of the con- 
ventional YM action. As it is symmetric, traceless 1 and 
divergence-free, it is the ideal candidate to describe a pure 
2 ++ state. But in our context we have mass scales, so the 
energy-momentum tensor, which differs from d9j, displays a 
nonvanishing trace, giving rise to a scalar state, so that we 
do not end up with a pure 2 ++ state. We could maintain |9]) 
as a candidate, but it is not conserved and so its divergence 
can create vector states. It is easy to construct a traceless and 
divergence-free rank 2 tensor from the operator F aiJ Fav- In 
particular, we propose here 



0, 



/'V 



with P pV = 5^ v - 



+d 2 ( 5„ v - \Pp> ) 3a3p9a|3 , 



3 2 



the transverse projector. By construc- 



tion, we have 0, 



0. 



■Vfl- 



next to d/j&fjv 



0. We also find 



®uu = 0, while for y = 0, ©^v = 3 4 A , V , so that 0^ v reduces to 
the derivative of the energy-momentum tensor. Thus, A , V ap- 
pears to be a useful local generalization of the normally stud- 
ied energy-momentum tensor. In addition, ©^v is identical to 
the rank 2 tensor employed in |fl9l . 

It remains to evaluate the spectral densities of the physical 
part of the scalar, pseudoscalar and tensorial correlation func- 
tions. For the threshold, we have x = 2(/j 2 + vV + 29 4 ) = 
1.3568 GeV 2 = With some calculational effort one finds 
(see also l|20l ) 



Jo V s 
with € {0 ++ ,(r+,2 ++ } and 

O 0++( 5 ) = c o++ 



E o / 1 4u 2 „ dt 

s v z — s 



{1+sT) 3 \2s 



29 4 



°o-+( 5 ) = c o-+ 



(1+sTf 

-7 



26 



4 , 



s 4s 2 



e 8 eV e 4 / % 2 3 



where the c^'s are positive constants, irrelevant for our pur- 
poses. The [0, 1] Pade approximants will be given by 



Po(z) 



-Vri 



(10) 



The trace anomaly is often ignored in a first order approximation (T). 
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Employing ( fT0| i, we immediately find a pole at z — v\ /v' , 
meaning that the mass estimate itself is given by mo++ = 
yVo/vp l n tne following figure, we show the results as func- 
tions of the subtraction scale T 

m(GeV) 




We now invoke the principal of minimal sensitivity to fix the 
scale T. Although we are unable to find a 7"* which would 
be a minimum or inflection point for m^(T), the situation 
changes when looking at the relative masses. We find an 
inflection point 2 in "' 2++ , resp. '" 2++ , given by T* m 0.34, 
resp. w 0.35. It is reassuring that this optimal T„ is al- 
most identical, while it is also located in the window where- 
fore niQ ++ < OT2++ < m Q-+- Setting T t . = 0.34, we come to 
main result of this letter, namely 

m 0++ « 1.96 GeV ,m _+ « 2.19 GeV ,m 2++ « 2.04 GeV . 

(11) 

Evidently, we should not have hoped to recover the lattice data 
very precisely, after all we are relying on a tree level propaga- 
tor. Though, it is interesting that we are able to reproduce the 
hierarchy mo ++ < nt2++ < mo h an d that the mass estimates 

2 We tried the same for "'°~ + , giving a zero in the third derivative, at T„ as 
0.35. 



( p"Tj ) are not too far from the lattice values HI 

moj + ~ 1 - 73 GeV » m o-+ ~ 2.59 GeV , mf ++ 2.40 GeV . 

Comparing our results to the lattice, one sees that we are 
within a 20% range. It is worth mentioning that we omitted 
instanton contributions, known to be relevant for the scalar 
and pseudoscalar channel, giving an attractive, resp. repulsive 
contribution around 200-300 MeV, see e.g. some works of JT]. 
We hope to come back to this issue in future work. 
Needless to say, many aspects of the Gribov-Zwanziger ap- 
proach are still open. For instance, the issues of how to give 
a precise definition of the physical subspace and how to prove 
its unitarity remain to be unraveled. Though, given the reason- 
able estimates for the masses we have obtained in a first sim- 
ple approximation, the study of the spectrum of the glueballs 
within the Gribov-Zwanziger framework looks very promis- 
ing and viable. Let us also mention that the IR moment prob- 
lem settled here could also be helpful to investigate the spec- 
trum of the glueballs in other approaches, also predicting a 
confining gluon propagator fl8l . One could try to get as good 
information as possible on the low momentum glueball opera- 
tors and feed this as input to the moment problem. Arguments 
based on IR power counting have already found applications 
in Schwinger-Dyson inspired studies, see e.g. |2D . We hope 
this letter will stimulate further research on glueballs, what- 
ever approach one is inclined to work with. 
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